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ABSTRACT 

We present the results of a pair of 100 ksec Chandra observations in the Small Magellanic Cloud 
to survey HMXBs, stars and LMXBs/CVs down to = 10^^ erg/s. The two SMC deep- fields are 
located in the most active star forming region of the bar. Deep Field-1 is positioned at the most pulsar- 
rich location (identified from previous surveys). Two new pulsars were discovered in outburst: CXO 
J004929.7-731058 (P=894s), CXO J005252. 2-721715 (P=326s), and 14 candidate quiescent pulsars 
were identified from their timing and spectral properties. Out of 12 previously known pulsars in the 
fields, 9 were detected, with pulsations seen in five of them. This demonstrates for the first time that 
a significant fraction (at least 60%) of these systems have appreciable accretion driven X-ray emission 
during quiescence. Two known pulsars in the field were not detected, with an upper limit of < 
5 xlO'^^. The full catalog of 394 point-sources is presented along with detailed analyses of timing 
and spectral properties. Future papers will report associated observations obtained with HST and 
Magellan to identify optical counterparts. 
Subject headings: Pulsars, X-rays 



1. INTRODUCTION 

The Small Magellanic Cloud has proven to be an in- 
comparable laboratory for astrophysics. The SMC is 
dwarf irregular satellite of the Milk Way, and unlike all 
other members of that group is experiencing an era of 
intense star formation. Situated in a part of the sky 
unobstructed by the galactic plane, the SMC affords a 
ringside seat to the unfolding drama. The combination 
of low extinction {nH ~10^^) and a small physical di- 
mensions in relation to its distance from earth (53 kpc ) 
effectively puts the entire population at a common dis- 
tance. A substantial fraction of the SMC can be observed 
simultaneously thanks to its compact size, which facili- 
tates population studies. 

Monitoring surv eys with e.g. RXTE 
(jLavcock et al. 20^51 ) only see pulsars in outburst 
at >10^^ yet most X-ray binary pulsars spend the 
majority of the time in quiescence. In order to obtain 
a measure of the true population, it is necessary to 
observe and count X-ray pulsars in their low luminosity 
state. The very concept of a quiescent low-luminosity 
state is an open hypothesis, so we target a densely 
populated environment at known distance and low 
nH with sufficient sensitivity (lO'^^erg/s) to detect the 
unseen majority. Based on observational evidence and 
theoretical calculations we predicted that Chandra could 
detect these systems, or place strong restrictions on the 
quiescent emission scenarios if none were found. 

The focus of this paper is the search for X-ray pul- 
sars in a quiescent state. That is, exhibiting emission 
(whether pulsed or not), at times other than during a 
type I or "normal" periastron outburst (and excluding 
type II, "giant" outbursts). In this work [5] we present 
new evidence for quiescent emission in previously known 



pulsars, whose expected outburst dates are predicted by 
ephemerides. However if the SMC harbors a large pop- 
ulation of systems that spend most of the time in quies- 
cence it will not be possible to find their orbital param- 
eters by X-ray monitoring. Instead these systems can 
only be found by deep pointed X-ray observations with 
sufficient angular resolution to guide searches for opti- 
cal counterparts. The orbital parameters can then be 
measured by photometric monitoring or radial velocity 
studies. There are two compelling scenarios that lead us 
to expect that such a population exists. 

(1) The physics of magnetic ally channelled accretion 
(see for example (|Frank et al.h ) predicts the existence 
of a centrifugal barrier. Magnetic field lines emanat- 
ing from the neutron star's polar regions co-rotate with 
the NS itself, causing plasma attached to the field lines 
beyond the co-rotation radius (where the keplerian or- 
bital period equals the NS rotation period) to be fiung 
away by centrifugal force. For accretion to occur, the 
ram-pressure of the infalling matter must drive the mag- 
netosphere boundary inside the co-rotation radius (rc). 
Matter orbiting inside Tc has angular velocity equal to or 
greater than co-rotation, and is channelled onto the NS 
polar caps, spinning up the pulsar in the process. This 
condition is a strong function of pulse period because the 
Tc gets progressively closer to the NS surface for increased 
spin rate. The most rapidly rotating neutron stars are 
therefor only expected to switch on as pulsars at very 
high mass transfer rates. Observational evidence for this 
picture is mounting, as many more outbursts are seen 
in SMC pulsars with periods longer than a few tens of 
seconds, compared to shorter period systems. 

(2) The majority of X-ray pulsars, and all but one of 
those in the SMC with identified optical counterparts. 
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have been fou nd in binary system s with Be-star com- 
panions (See (jMcBride et al. 2008^ for a recent review 
of the evidence). The accreted material powering X- 
ray emission from the pulsar comes from a circumstel- 
lar disk around the Be star. The exact mechanism 
of disk formation is unclear at present, but progress 
has been made in understanding the disk-NS interaction 
([Okazaki and Negueruela. 20011 ). Several hues of obser- 
vational and theoretical evidence point to the role of or- 
bital eccentricity in producing regular "normal" (type-I) 
X-ray outbursts in Be-HMXB. Systems with circular or- 
bits (or nearly so) are thought to exhibit only occasional 
giant (type-II) outbursts. It is further speculated that 
these systems can experience constant low-level accre- 
tion, yielding X-ray luminosites not detectable by the 
pre-Chandra generation of instruments. In the Milky- 
Way, the Be/X system X Persei is the prototype of this 
class with eccentricity=0 .11, and there are at least two 
other galactic examples (|Delgado-Marti et al. 200"ll) . X 
Per fortuitously lies somewhat out of the galactic plane 
otherwise it would be hidden by dust/gas, as we presume 
are many of its brethren. The X-ray luminosity of X Per 
varies, but is seldom above Lx=lO^^ erg/s. 

Our expectation is that X-ray binaries containing fast- 
rotating neutron stars, and those with circular orbits are 
under-represented in current catalogs. In both cases due 
to low quiescent luminosity and the infrequent occur- 
rence of giant outbursts. One of our aims is to provide a 
complete sample of X-ray binaries, with greatly reduced 
luminosity bias compared to the existing sample. Stel- 
lar population synthesis models will greatly benefit from 
hard numbers upon which to base the relative probabili- 
ties of various evolutionary paths. For example in mod- 
elling the supernovae that occur in formation of XRBs, 
one needs to reproduce the actual distribution of orbital 
periods and eccentricities. If current catalogs under rep- 
resent systems with large orbital separations and low rj 
the input to these models will be flawed. 

Studies of a rchival Chandra da ta for a large number 
of galaxies by (|Grimm et al. 20031 ) have found an appar- 
ently universal correlation between star formation rate 
(SFR) and the integrated luminosity of HMXB. The 
absolute number normalization of this relationship is 
however unknown, because sufficiently faint fluxes have 
not been probed to discover the underlying numbers of 
HMXBs. It is for example physically relevant to know 
the number of neutron stars formed in binaries per unit 
SFR, which one might not be able to tell purely from the 
luminosity function, if it is missing a significant sub-class 
of HMXB. 

We present the complete X-ray source catalog as an 
electronic table accompanying this paper. We describe 
below in detail the properties of the brightest sources, 
and give an analysis of the timing and spectral properties 
of the sample, focussing on X-ray p ulsars. An optical sur- 
vey is being conducted in parallel (jAntinou et al. 20081 ) 
to enable more robust and complete classifications of X- 
ray binaries. Cataclysmic variables, stellar coronae and 
background acti ve galactic nuclei (AGN). Initial results 
are reported by (jZezas et al. 2008t ). 

2. OBSERVATIONS AND DATA REDUCTION 

The two SMC deep-fields are located in the most active 
star forming region of the bar. Deep Field-1 (DFl) is po- 



sitioned at the most pulsar-rich location identified from 
previous surveys. The DFl coordinates also coincide 
with the primary position used by the RXTE monitoring 
program (jLavcock et al. 2005f) (herafter L05), for weekly 
pulsar observations of 5-10 ksec duration for the past 8 
years. Deep Field-2 (DF2) was previ ously observed wit h 
Chandra in 2002, as reported by (|Tavlor et al. 2003f) . 
The earlier observation was a 10 ksec ACIS-I observation 
which discovered two new pulsars and a bright supernova 
remnant. The two SMC Deep Fields were each observed 
for a total of 100 ksec as detailed in Table [TJ For mission 
operational reasons the observations could not be sched- 
uled as contiguous blocks, a circumstance anticipated in 
advance. DFl was observed in two observations of ~ 50 
ksec each, over the course of two days. DF2 was split 
into 3 observations spread over three days. Aim-point 
and roll-angle were identical for the observations mak- 
ing up each field, so as to simplify the data reduction. 
The Advanced CCD Imaging Spectrometer (ACIS) was 
used, with the ACIS-I array at the aimpoint and CCDs 
10-13 & S5-S6 taking data in the full unfiltered sensitiv- 
ity range (nominally 0.3-8 keV). Timed exposure, faint- 
mode was used with the default "most-efficient" frame 
time of 3.2sec. Smoothed images of the Chandra data 
are shown in Figures [T] & [2l 

The data were stacked and reduced using the process- 
ing pipeline developed at CfA for the Chandra Multi- 
wavelength Project (ChaMP, (|Kim et al. 20041) ). The 
pipeline performs all standard Chandra calibrations and 
corrections including exposure-map correction, cosmic- 
ray afterglow removal, aspect solution, charge-transfer 
inefficiency, effective area, and time-dependent energy re- 
sponse (due to contaminant buildup on the ACIS optical 
blocking filter). Source detection was performed with 
the wavedetect algorithm, on images extracted in three 
energy bands (Bx=0.3-8 keV, Sx=0.3-2.5 kcV ,Hx^2.5- 
8 keV) at all available spatial scales (1,2,4,8,16). X- 
ray source properties were measured using the apertur e 
photometry procedure described by (|Hong et al. 20051 ). 
which accounts for neighboring sources with overlapping 
source or background regions. The full position depen- 
dent instrument sensitivity is accounted for by applica- 
tion of the exposure-maps and response matrices gen- 
erated by the ChaMP pipeline. Positional error circles 
are assigned to each source based on net counts and dis- 
tance fro m the aim-point, ac cording to the now standard 
model of (Hong et al. 2005h which is based on detailed 
simulations of the full Chandra optical path and ACIS 
characteristics. Source detection and photometry were 
performed both on the observations in table [T] and the 
complete stacked data-set for each field. The resulting 
object lists were then combined to yield a catalog of 
unique sources, with the best measurement of each source 
used to define its coordinates and error circle. In nearly 
all cases the stacked Bx image provided the best mea- 
surement with the following exceptions. Variable sources 
which may be detected at better S /N in a particular ob- 
servation, and particularly hard or soft sources which 
may achieve the best S/N in the Sx or Hx images. 

3. TIMING ANALYSIS 

Light-curves were extracted in units of net rate (i.e 
background corrected counts per second) using good 
events from source and background regions of every suf- 
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ficiently bright object detected. The hghtcurves were ex- 
tracted following the usual CIAO analysis thread with 
the tool dmextract at time resolutions of 1, 4, 8, & 
16 times the intrinsic ACIS frame-time of 3.2s. Power 
density spectra (PDS) were computed using the Lomb- 
Scargle method ( Scargi e 1982) and all periods above 90% 
significance were recorded. Due to scheduling constraints 
our 100 ksec observations were split (see table [T]), so the 
lightcurves were analysed both independently and whole. 
The number of counts needed to unambiguously detect a 
pulse period depends on the ratio of period to observation 
duration, sampling interval, pulse fraction and pulse pro- 
file. Taking an optimistic cut at >100 net counts in 100 
ksec, provides a lower threshold of approximately 10*^4 
erg/s for detecting periodicities in SMC pulsars. The 
threshold for detecting a source without regard to timing 
properties is about ten times fainter. For comparison the 
pulsation detection thre shold for the RXTE SMC pul- 
sar monitoring program (|Lavcock et al. 20051 ) is approx. 
IQ^s erg/s, and only pulsed sources can be detected. A 
periodogram analysis was conducted blind (period search 
range 6-2000s) on every source with more than 50 net 
counts per lightcurve segment, or 100 net counts for the 
stacked data. Deep Fields 1 and 2 contained 19 and 
14 sources respectively with >100 net counts. Spurious 
periods are expected near the Chandra dither periods 
(1000s, 707s), their harmonics (principally 500s, 353s) 
and the ACIS sampling period (3.2s). Timing artifacts 
are only expected for sources that fall close to a CCD 
edge or node boundary. In this situation the dither off- 
sets can periodically move all or part of the source region 
on and off the detector, resulting in a modulation of the 
count-rate. Examples of dither artifact can be seen in 
the PDS of several sources and are identified in the cap- 
tions and/or accompanying text. For each source with 
a positive period-detection, the data were folded at the 
most probable period and grouped inlO bins per period. 
In order to analyze the pulse profiles in a uniform man- 
ner, our folding code automatically establishes the epoch 
of peak flux by folding the data using the first data-point 
as a trial epoch. The data are then folded a second time 
using the derived epoch. Pulsed flux fractions were com- 
puted from the binned pulse-profiles according to the for- 
mula fp = {fmax - !min)l !max- We notc that /p teuds 
to unity as the un-pulsed fiux approaches the sensitiv- 
ity limit of the observation. Thus pulse fractions have a 
natural tendency to be overestimated and values consis- 
tent with 1.0 are in fact upper limits. The phase- folded 
data-points were also process ed with the loes s algorithm 
(a non parametric smoother () Cleveland 19811 )) to obtain 
an un-binned visualization of the pulse profile. 

4. NEW PULSARS 

Two previously unknown pulsars were discovered in 
outburst with periods of 892s, 326s, and fourteen quies- 
cent pulsar candidates were identified, principally from 
their timing properties. 

4.0.1. CXO J004929.7-731058 - SXP892 

Detected in Deep Field 2 with 1434 net counts. Fig- 
ure m the 892s pulsation is detected at very high signif- 
icance as shown in figure [3l The derived luminosity is 
Lx = 9.3x10'^^ erg/s, consistent with a faint type-I (nor- 
mal) outburst from a Be X-ray binary pulsar. The ACIS 



spectrum from 0.35-8 keV is well fit by a power-law with 
photon index 7=l[3l with very low extinction. DF2 was 
observed by Chandra in 2002 for 10 ksec, SXP326 was not 
detected. Had it been active at the same level it would 
have registered ~143 counts, thus it's transient nature is 
confirmed. It has also never been seen by RXTE, which 
has occasionally observed at positions with significant 
response at the coordinates of DFl. 

4.0.2. CXO J005252.2-721715 - SXP326 

Found in Deep Field 1 with 9292 net counts. The 326s 
pulsation is detected at very high significance as shown 
in figure m The derived luminosity is = 5.5x10'^^ 
erg/s, consistent with a moderate type-I (normal) out- 
burst from a Be X-ray binary pulsar. The ACIS spectrum 
from 0.35-8 keV is well fit by a power-law with photon 
index 7=l[3l with very low extinction. DFl is monitored 
regularly by RXTE, which has never detected SXP326, 
therefore either the pulsar is rarely in outburst, or its 
normal outbursts are always below the RXTE detection 
limit of ~10^^ erg/s 

4.1. Candidate Quiescent Pulsars 

We show in Table [2] the timing- analysis parameters for 
14 faint pulsars discovered in this work. The luminosities 
of these sources are in the range of a few 10'^3~34 erg/s, 
substantially below the level of a normal (Type-I) out- 
burst (>10'^^ erg/s). It is also unlikely that so many sys- 
tems would be at periastron simultaneously. We there- 
fore propose we are detecting pulsars with persistent low 
luminosity X-ray emission. The pulse periods of these 
sources (6.3s - 655s) span the usual range. Comparison 
with the period distribution of pulsars seen by RXTE 
shows a slight excess of quiescent pulsars with periods 
below 100 sec, but a KS test shows it is not significant 
(P=0.34). All but one pulsar were found in the ACIS- 
I data, the source XS07156B7_001 on the ACIS-S array 
was found with a period 6.3s. 

Two supernova remnants (SNR) are visible in Fig- 
ure [21 and each has a candidate pulsar within its 
borders. CXO J005059.0-732055 with P=9.39s is in- 
volved in th e outer region of the bright SNR dis- 
covered by (iTavlor et al. 2003D . The 8.94s source 
CXOPSJ004905.1-731411 (Figure[T2]) appears to be near 
the center of the fainter SNR. These alignments could be 
coincidence given the density of sources in the field and 
size of the SNRs. 

The source CXO J005215.4-731915 (Figure ^ was 
detected in Deep Field 2, and lies within 7.5 arcsec of 
SXP15.3, however we detected a seperate and distinct 
source at the exact reported coords of SXP15.3, as can 
be seen in Figure [2l What appears to be the same 
object was detected in one of the individual exposures, 
as XS08479B2_014 with pulse period 8.234s at 94.2% 
significance so there is some ambiguity as to the correct 
pulse period. 

The spectral/quantile analysis presented in section [7] 
identified three sources with enough counts for unam- 
biguous determination of their spectral index. In the 
quantile plane (S/N/^10) they are indistinguishable from 
known pulsars. Power spectra of these sources were 
re-examined for weak pulsations. For CXO J005057.1- 
731008 (Figure [9]) we found a period of 98s (71% sig- 
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nificance). For CXO J005215.4-731915 (Figure we 
found a weak periodicity at 7s (62 % significance). In 
CXO J004905.I-73I4II (FigurelHl) we found a period of 
8.94s (93.3% significance). A lower detection threshold 
is justified because there is a strong reason to suspect the 
sources are pulsars based on the spectral evidence. The 
detection threshold used for our other pulsar candidates 
was set at 90% as a result of these spectrally supported 
identifications. 

5. KNOWN PULSARS 

There are 12 known X-ray pulsars covered by these 
observations: seven in Field-1 and five in Field-2. We 
detected strong pulsations in six of the known pulsars 
(SXP756, SXP316, SXP172, SXP138, SXP59, SXP8.89) 
and tentatively detect weak periodicity in one more 
(SXP15.3). Two additional pulsars (SXP9.12, SXP7.78) 
were detected but without pulsations being discernible. 
Most of the known Be/X-ray binaries in the SMC 
are seen regularly b y the RXTE m onitoring program 
(|Lavcock et al. 2005f ). (jGalache et al. 2008) . and as a re- 
sult have known X-ray ephemerides. We calculated the 
X-ray outburst phase for each detected pulsar in order to 
search for quiescent emission in systems far from perias- 
tron. The ephemeris convention used by LOS places the 
peak in the folded X-ray lightcurve at phase=0, while 
G08 place the peak at phase=0.5. The RXTE data pro- 
vide additional because we know the general state of ac- 
tivity in each pulsar around the times of the Deep Field 
observations. 

Three known pulsars in our fields were not detected: 
SXP4.78 has no ephemeris. SXP46.6 was at phase 0.33 
(G08 ephemeris, P=137d) which is sHghtly before flux 
is typically detected by RXTE . SXP82.4 was at phase 
0.86 (G08 ephemeris, P=362d) which is after the normal 
phase-range of RXTE detections. 

5.0.1. CXO JOO4942.O-732314 - SXP755 

Detected with pulse period 746s in deep field 2, the 
outburst phase (LOS ephemeris) of the system was 0.7, 
in the quiescent portion of the long-term lightcurve, and 
presumably far from periastron. 

5.0.2. CXO J005044.6-731605 ~ SXP323 

We detected a bright source with pulse period 316s, 
located at the coordinates of SXP323, it is likely 
the same pulsar. Although the period differs signifi- 
cantly from the long-time average value of 323s, the 
RXTE monitoring program recently detected the sys- 
tem in outburst on MJD 52960 with pulse period 317s 

100% period-detect ion significance, uncertainty<0.1s). 

Galache et al. 20081) attributed that outburst to a Type 
II (super-outburst) event, which we suggest has spun-up 
the pulsar. At the time of the Deep-Field-2 observation 
(MJD 54061) the outburst phase of SXP323 was 0.8 (G08 
ephemeris) and the pulsar was far from periastron. 

5.0.3. CXO J005151. 9-731033 - SXP172 

SXP172 was detected with pulse period 171.8s in Deep 
Field 2 (MJD 54061 ) with 2-10 keV luminosity L^ = 
l.S3x lOSSerg/ s. The pulsar has historically shown a pro- 
longed episode (MJD 51S00-S2S00) of steady spin-up (see 
(jGalache et al. 2008f) . figure 30). The period detected 



here is in line with the end of that active episode, in- 
dicating the pulsar has not continued to spin up during 
the ~1000 days since its activity dropped below the de- 
tection threshold of RXTE, nor has it spun down appre- 
ciably The X-ray outburst phase for SXP172 on MJD 
54061 was 0.32 (LOS ephemeris), placing it reasonably 
far from peak-phase. This pulsar has historically been 
seen in outburst at times uncorrelated with the X-ray 
and optical period. 

5.0.4. CXO J005323.9-722715 ~ SXP138 

This source is confirmed by its coordinates as SXP138, 
in addition to the previously known pulse-period of 
138s we detected a number of other periods. Obser- 
vation XS071SSB2: 688 net counts. Periods: 972.76s, 
500. 2Ss, 333.33s , 138.88s (972/7), 83.97s, Observation 
XS07327B2: 537 net counts. Periods: 970.87s, 334.45s, 
161.29s (970/6), 138.89s (970/7). We note the 138s pe- 
riod contained the most power as shown in figure[8l Anal- 
ysis of the combined light-curve of both observations 
revealed highly significant periods at 500.000, 333.333, 
161.290, 138.888. The exact factors (500, 333.33) of the 
Chandra Y-axis dither period (1000s) suggests at least 
some of these periods are due to the source being moved 
back and forth across a chip edge or node-boundary. In- 
deed CXO J005323.9-7227151ies within 10" of the CCD's 
edge. It is not easy to see why the remaining periods 
(972,161,138) should follow a separate precise relation, 
unless one of them is a pulse period and the others are 
harmonics. Careful analysis of the two individual ob- 
servations and the full lightcurve shows that although 
both segments alone show a 138s period modulated at 
972s, the phase of the long periodicity is different in the 
two obsids and in consequence is totally suppressed in 
the full lightcurve. We conclude the additional periods 
are an interaction between the pulse period and Chan- 
dra's dither. At the time of the Deep-Field-1 observation 
(MJD 53851) SXP138 was at outburst phase 0.005 (G08 
ephemeris), which is the point of minimum flux, and as- 
sumed to be far from periastron. 

5.0.5. CXO .1005456.1-722648 - SXP59 

The source was detected in the Chandra obser- 
vation of Field 1 at a luminosity of 1.83x10'^'' 
erg/s with p ulse period 58.8 s. According to the 
ephemeris of (i Galache et al. 2008D . MJD 52306.13.7 -\- 
nx 122.10±0.38 d, which defines phase=0 to coincide 
with peak X-ray luminosity in the lightcurve folded at 
the orbital period. SXP59 was at outburst phase 0.65 
during the Deep Field 1 observation on MJD 53851. We 
measured a pulse period of 58.85815s whi ch is longer than 
the last reported measurement by f Gala che et al. 2008f) 
on MJD 54000, but in line with the long term aver- 
age value. SXP59 tends to spin-up when in outburst, 
and spin-down in quiescence. The outburst phase (G08 
ephemeris, P=122d) was 0.6 during the DF-1 observa- 
tion, which is within the range of normal outbursts seen 
in SXP59. 

5.0.6. CXO .7005213.9-731918 - SXP15.3 

Deep Field 2 registered a clear detection at the ex- 
act coordinates of SXP15.3. We detect 147 net counts 
{Lx — 6.77 X lO^'^ erg/s) with a powerlaw spectrum. We 
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do not see any significant p ulsations. T here is no X-ray 
ephemeris for SXP15.3, but (|Edge 20051 ) produced an op- 
tical ephemeris (P=75d) from optical photometry of the 
counterpart. Accordingly the (optical) outburst phase 
was 0.13 on MJD 54061. 

5.0.7. CXO J004913. 5-731138 - SXP9.13 

We detected a bright (1178 net count) source in DF- 
2 at essentially identical coords to those reported for 
SXP9.13 (within our 0.5' error radius). However pulsa- 
tions were not detected, which is remarkable given the 
source brightness. The time resolution of the Chan- 
dra ACIS data is 3.2 sec, which is sufficient to criti- 
cally sample the 9.13s pulse period, the pulse fraction 
must therefore be very small. We measured a luminosity 
of \jx = 8.5 X lO^"* erg/s (from fluxBc). The outburst 
phase (G08 ephemeris, P=77d) was 0.82, which is out- 
side the range of the normal outbursts seen in SXP9.13, 
and therefore presumably far from periastron. 

5.0.8. CXO J005153.3-723148 - SXP8.88 

Detected in Deep-Field 1, with pulse period 8.89s. The 
outburst phase was 0.75 according to the ephemeris of 
([Corbet et al. 20041 ) (P=28d) who defined peak flux at 
phase=0. This detection is therefore outside of the nor- 
mal range of outbursts. The source fell on the ACIS-S 
array observation of DFl (XS07327B2_024). 

5.0.9. SXP7.78 

A source was detected in Deep field 1 at the coordinates 
of SXP7.78. With 72 net counts the detection is highly 
significant and yields accurate position and quantile in- 
formation, but the lightcurve is insufficient for period 
analysis. The outburst phase (G08 ephemeris, P=43d) 
was 0.62 which is within the normal range for SXP7.78 

6. VARIABILITY OF NON-PULSED SOURCES. 

X-ray variability is the hallmark of accretion and is 
also characteristic of fiaring stellar coronae. Cataclysmic 
variables (accreting WDs) and Low-Mass X-ray binaries 
(LMXBs: both NS and BH systems) exhibit aperiodic 
variability (flickering) over a wide range of timescales 
from milliseconds up. Such variability produces a char- 
acteristic power-law distribution (red noise) in the PDS 
which can be spotted by visual inspection. 

PDS analysis is not feasible with fewer than about 100 
counts, so Ughtcurves binned by 16 times the ACIS frame 
time were visually inspected for all sources with at least 
50 counts in the stacked dataset. This search turned up 
a large flare in CXO J005428. 9-723107, lasting for -3000 
sec at the end of the first of two 50 ksec exposures of Deep 
Field 1, shown in figure [TBI None of the other 95 sources 
so examined showed any similar feature. The power spec- 
trum of the flare is shown figure [T7[ no periods, QPOs, 
or red noise are evident. Subtracting the total net counts 
for the source during the second exposure from the first 
(259cts - 32cts), we obtain an estimate of 277 counts in 
the flare. Comparing quantiles for the source between 
exposures, the flare (median energy E50 = 1.288 keV) 
was harder than the quiescent emission (E50 = 1.021 
keV). The quantile values during and post flare are plot- 
ted in the bottom panel of Figure 1141 During quiescence 
the source is consistent with a fairly cool thermal plasma 



(kT'--^ 0.4 keV) and the extinction is clearly low although 
poorly constrained. During the flare the quantile moves 
to higher temperature (kT~lkeV) and we obtain a better 
constraint on the extinction, confirming its very low value 
(nH<10^"'^ cm). A fit to the X-ray spectrum of the ob- 
servation containing the flare is shown in flgure 1181 The 
fitted model is a Raymond-Smith model (optically thin 
thermal plasma) with solar abundance, plasma tempera- 
ture kT=2.5±0.4 keV and essentially zero extinction, the 
reduced — 1-14 (20 degrees of freedom). Fixing the 
extinction at nH=10^^ makes the slightly worse (1.4). 
The Mekal model gave consistent results (kT=1.7±0.3 
keV, x^—1.2). For comparison we fit an absorbed power- 
law (r=3±0.3, nH=2.6(±l)xl02i x^=0.64). 

7. SPECTRAL ANALYSIS 

X-ray spectral analysis requires a reasonable number 
of photons, and was performed for the 14 sources with 
more than 250 net counts in the stacked dataset. Spectra 
were extracted from the XPIPE source and background 
event files using CIAO 3.4, and grouped to give 10 source 
counts per spectral bin. Spectral fits were performed 
with XSPEC, the results are given in table [3l All but 
one of these sources are pulsars and are well described 
by the power-law model typical of X-ray pulsars over this 
energy range. Spectral indices are 0.7-1.8 and derived 
neutral hydrogen column densities lie in the range 0.6 
- 10 xlO^^. Pile-up (multiple photons arriving within 
a single ACIS pixel during one 3.2s readout.) can be 
safely ignored for all our sources, the brightest of which 
contributes 0.10 counts/sec 

The majority of the sources in our catalog are too 
faint for spec tral fitting, we t herefore turn to quan- 
tile analysis (|Hong et al. 2004 ). Quantile analysis is a 
method to optimally extract spectral information from 
low-count X-ray sources, without imposing arbitrary 
energy-bands choices. The advantage over X-ray colors is 
that all sources can be represented in the same parame- 
ter space, with a somewhat uniform error distribution. 
Energy band choices inevitably lead to some sources 
having a gross imbalance in the nu mber of counts in 
each band. Following the methods of (jHong et al. 2005f) . 
background corrected quantile QDx and QDy values were 
computed from each source's photon distribution (event 
list), and spectral model grids were generated using the 
full Chandra ACIS-I response function. The quantiles 
values with grids for absorbed power-law and thermal 
bremsstrahlung are plotted in figures fT4l fc [T5l for a series 
of cuts in signal-to-noise-ratio (S/N> 3, >5, >10). In 
both fields, all of the confirmed pulsars (both with and 
without detected pulsation in the Chandra data) clus- 
ter around power-lawr=l, and neutral hydrogen column 
density nH <10^^. Several unidentified sources lie in this 
region also, and are therefore consistent with being addi- 
tional X-ray binary pulsars. In particular there are two 
objects with S/N>10 (one in Field 1 & one in Field 2) 
whose quantiles are indistinguishable from specific con- 
firmed pulsars within the errors. We therefore suggest 
the following objects are pulsars: CXO J005057. 1-731008 
(XS58479B3-003, net counts = 700), & CXO J005215.4- 
731915 (XS58479B3-022, net counts= 286). These four 
sources were further scrutinized for timing features below 
the 90% blind search significance described in section [31 
This analysis revealed pulsations in the two brightest 
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sources CXO J005057.1-731008 & J005215.4-731915 (see 
section 11]). A group of sources with very large extinc- 
tion occupy the upper right corner of the power-law grid. 
These are likely obscured active galactic nuclei (AGN) 
with significant internal extinction. The quantile pa- 
rameters for these AGN-candidates are consistent with 
r ^^1.7 power-law spectra, typical of the class. 

The upper left region of the quantile diagram contains 
a group of soft sources that extend outside the power- 
law grid, and are consistent with thermal models ap- 
propriate for stellar coronae. In the figures we show a 
model-grid for thermal bremsstralung, and the sources 
occupy a region approximately bounded by plasma tem- 
peratures in the range kT = 0.4 ~ 2 keV and nH <10^^cm. 
We note that grids can be constructed for more sophis- 
ticated stellar-coronal models such as Raymond-Smith 
or MEKAL, but for low count sources they provide no 
additional information. 

After accounting for stars, obscured AGN, and con- 
firmed pulsars, we are still left with the majority of our 
sources. The unknown sources are not randomly dis- 
tributed in the quantile-plane but lie predominantly in 
the region bounded by nH<10^^ and 2 < F >1 suggesting 
they are mostly accretion powered sources in the SMC 
(X-ray binaries, CVs) and beyond (AGN). X-ray binary 
pulsars can exhibit relatively soft quantiles with F ~ 2, 
see for example figure 3 in (jMcGowan et al. 2007| ). so 
it is plausible for some of these objects to be pulsars 
observed in a low luminosity or quiescent state. Cat- 
aclysmic Variables (CVs) and Low-Mass X-ray binaries 
(LMXBs) are also expected to appear in this region of 
the quantile diagram. Additionally CVs and LMXBs are 
likely to be faint with 0.5-10 keV luminosities in the range 
Iq31-35^ which is consistent with many of these sources. 
However, the greater proportion of the objects in this 
group are likely to be background AGN, which are typ- 
ically characterized by power-law spectra with F ^--^1.7 
and low intrinsic absorption. The expected nH distribu- 
tion for background AGN follows the line of sight column 
through the SMC. 

Determination of the nature of each individual source 
is being sought by an optical counterpart survey 
(|Antinou et al. 20"08[ l. based on photometry obtained 
with HST and the 6.5m Magellan-Baade telescope. 

8. CONCLUSIONS 

The complete source catalog is available on-line in elec- 
tronic form, it contains positions, count-rates, fluxes and 
quantiles for 394 X-ray sources (211 in DFl & 183 in 
DF2). 

We have detected 16 X-ray pulsar candidates in the 
SMC Chandra Deep Fields (7 in DF-1, 9 in DF-2) based 
primarily on periodogram analysis of the lightcurves. 
Two of the newly discovered pulsars were in outburst, 
and their periods were measurable to a high degree 
of accuracy. We have therefore assigned the names 
SXP892 & SXP326 following the convention used by t he 
SMC RXTE monito ring project (jLavcock et al. 2005( 1 , 
(jCalache et al. 2008^ . An additional three candidates 
were detected with sufficient S/N to place them unam- 
biguously on the quantile diagram in the region inhabited 
by pulsars, in tight agreement with the powerlaw spec- 
tral model (F=l, nH^-lO^^) characteristic of SMC pul- 
sars. The remaining 11 pulsar candidates were identified 



from periodogram detections with blind search signifi- 
cance levels between 90%-95%, or in sources with fewer 
than 250 net counts. We report these candidate pulsars 
without an SXP name due to the uncertainty attached 
to their period determinations, for example we cannot 
be confident that the fundamental is the most significant 
peak in the periodogram, or that a small fraction (~10% 
for the 90% significance threshold) of these periods could 
be spurious. In the following discussion we assume that 
all the candidate pulsars are real detections. 

There were 12 X-ray pulsars known in the survey re- 
gion prior to these observations, 7 in DFl and 5 in DF2. 
Thus we have increased the sample of pulsars in the re- 
gion to a total of 28 {[7+7] + [5+9]). Most interesting 
is the fact that we were able to detect all the known 
pulsars in DF2 and most (4/7) in DFl. One of the pri- 
mary objectives of this project was to find out if HMXB 
pulsars could be detected between outbursts, and to de- 
termine their quiescent luminosity. Using the latest X- 
ray ephemerid es produced b y the RXTE monitoring pro- 
gram (Galac he et al. 20081 ) we can calculate the (pre- 
sumed orbital) phase of each pulsar at the time of the 
Chandra Deep Field observations. The ephemerides were 
constructed from pulsed-flux lighturves spanning about 
10 years of weekly RXTE observations. As such they 
predict the most likely time of X-ray outburst, which 
is presumed to reflect orbital modulation of the mass- 
transfer rate from the mass-donor star to the neutron 
star (the actual orbital ephemeris is not known for any 
SMC pulsars except SMC X-1). Of the 10 pulsars with 
an X-ray ephemeris fromn RXTE, 6 were detected far 
from their normal range of outburst-phase and with lumi- 
nosity well below outburst levels. Two pulsars (SXP59, 
SXP7.78) were detected during their normal range of out- 
burst phase. Three known pulsars were not detected, 
which for a requirement of 5 counts in lOOksec gives an 
upper limit of < 2.6x lO^^erg/s. Of the undetected 
pulsars, two have X-ray ephemerides and were both ob- 
served outside their normal range of outburst-phase. 

By detecting 6/10 pulsars for which there are 
ephemerides, outside of their predicted outburst phase, 
we have demonstrated that a significant fraction (at least 
60%) of these systems have appreciable accretion driven 
X-ray emission during quiescence. To the six confirmed 
quiescent pulsars we can add the 14 candidate quiescent 
pulsars, whose 2-10 keV luminosities are substantially 
lower than normal Be-pulsar outbursts. This new sam- 
ple is consistent with a large population -about twice the 
size of the known pulsar population- of quiescent "nor- 
mal" pulsars. It is also possible that some are in a per- 
petually low luminosity state, such systems are known in 
the Galaxy, for example X Per. The theoretical work of 
(jOkazaki and Negueruela. 20011 ) informs us that Be/NS 
systems with low orbital eccentricity accrete at a very 
low constant rate. The Be star's circumstellar disk in 
such systems is permanently tidally truncated (at the 3:1 
resonance) well within the neutron star's orbital radius, 
except during large mass-ejections (the assumed cause of 
type II "giant" outbursts). This picture neatly explains 
how more than 50% of the X-ray pulsar population has 
remained hidden from RXTE and is only now coming to 
fight. 

Pulse profiles created for every pulsar detected in the 
Chandra data show a range of forms. The typical pro- 
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file shows significant flux above the un-pulsed component 
for most of the pulsation cycle. A qualitative distinction 
exists between profiles that are broad and asymmetric 
(e.g. SXP326), and those that are strongly peaked which 
tend to be narrower and more symmetric (e.g. SXP172). 
These differences are likely due to differences in geome- 
try of the emission regions (polar caps). The literature 
documents several efforts to model X-ray pulse profiles in 
terms of combinations of fan and pencil beams (See for 
example (jParmar et al. 1989f )). Broad and cuspy pro- 
files are attributed to fan beam geometry, where the X- 
ray beam has the form of a (somewhat irregular) hollow 
cone, due to obscuration of the central part of the polar 
cap by the incoming accretion stream, or reflection/re- 
processing effects. Narrow profiles more closely resem- 
ble radio pulsars and are accordingly interpreted in the 
"light-house" model, as coUimated beams (pencil-beams) 
emanating directly from the polars caps. 

Including only sources with S/N>5, and counting the 
sources populating regions defined on the quantile dia- 
gram (Figures [Tib & [TSb) , we estimate of the propor- 
tion of candidate HMXB, Stars, AGN & obscured AGN 
in each field. These numbers are then used to estimate 
the absolute number of XRB, stars & AGN in the Deep 
Fields. We multiply by the fraction of high S/N sources 
in each category by the total size of the catalog down to 
the completeness limit of our dataset. Uncertainties are 
computed by taking ^/N of the number of 5ct sources, 
propagated to the full sample. We defined a flux limit 
corresponding 10 count source at the Chandra/ ACIS-I 
aimpoint in our 100 ksec stacked dataset, to approxi- 
mate the completeness limit of our sample. We do this 
calculation using a flux limited sample because our flux 
calculations include correction for Chandra's changes in 
sensitivity and spectral response with distance from the 
aim-point. If a count limit were used instead, it would 
not be uniform across the field. Each source has a 0-5-8 
keV flux computed based on its quantile values and the 
exposure map, a 10 count limit at the aimpoint translates 
to a flux hmit >10~^^ erg/cm^/s. 

In Deep-Field-1 we find 10.7% Stars, 53.3% AGN (of 
which 9% are obscured), and 35.7% XRB. Mufliplying 
these percentages by the number (131) of sources down 
to our completeness limit yields 14±6 stars, 71±17 AGN, 
and 46±13 HMXB. 

The same analysis of Deep-Field-2 yields 10.5% Stars, 
50% AGN (of which 7.9% are obscured), and 39.5% 
XRB. Multiplying these percentages by the number (107) 
of sources down to our completeness limit yields 11 ±6 
Stars, 54±16 AGN, and 42±14 XRB. 

For comparison with our observations, we first es- 
timated the expected number of bac kground AGN i n 
our survey based on the work of (iKim et al. 20071) . 
For the same fl ux limit of >10~"'^^ erg/cm^/s, 
(|Kim et al. 2007fl Fi gure 4 shows of order 10'^ AGN per 
sq. degree. Thus for an ACIS-I field of view (16'xl6') 
one expects to see ~71 AGN per field. The purpose of 
this comparison is to determine whether our proposed 
large expansion in the X-ray binary population of the 
SMC is reasonable in the light of what is known about the 
background AGN population. The AGN number for DFl 
is remarkably consistent with our prediction b ased on the 
AGN luminosity function of ()Kim et al. 20071 ) . We note 
that DF2 contains a large X-ray bright supernova rem- 



nant which likely explains the slightly (23%) lower AGN 
numbers in this field. Systematic errors were estimated 
by varying the boundary of the region in the quantile 
plane used to classify the AGN. Moving the boundary 
between r=1.2 -1.5 changed AGN vs XRB fractions by 
an insignificant amount compared to the statistical error. 

Based on the above results, we argue for a large qui- 
escent population of High Mass X-ray binaries (HMXB) 
in the SMC. Only a small fraction of these systems are 
in outburst at a gi ven instant in time. The work of 
(|Grimm et al. 2003[ ) demonstrates a universal relation 
between the integrated X-ray luminosity of HMXBs in a 
galaxy and the instantaneous star formation rate. This 
result now must be seen in the light of a similarly uni- 
versal duty cycle distribution for HMXB populations. 
At any instant, the proportion of systems in outburst 
is roughly the same. Although more numerous, the qui- 
escent systems evidently do not contribute significantly 
to the total luminosity of the galaxy. 

We note that Cataclysmic variables and LMXB are ex- 
pected to occupy the same quantile parameter-space as 
the AGN, due to their soft power-law and/or hot ther- 
mal spectra. Given that our AGN-count is consistent 
with predictions, there is not a great deal of room for 
additional source populations. Thus we can already con- 
strain the total number of LMXB and CVs to be within 
the error bars of the AGN estimate. 

The further goals of this project are as follows, and are 
being pursued in a series of papers in collaboration with 
our co-proposers. (1) Classify all the sources by find- 
ing and studying their optical counterparts. Discrimi- 
nating the different source classes is a requisite for aims 
3-5 below. Our optical survey is being carried out in 
parallel with the X-ray survey using the Hubble Space 
Telescope and Magellan 6.5m telescope. With high reso- 
lution imaging/photometry to R,I<26 we will be able to 
independently classify the optical counterparts to most 
of the sample. The optical counterpart identifications 
will enable items 2-4 below. 

(2) Construct the luminosity distribution down to faint 
X-ray fluxes, this will test whether multiple point source 
populations may be present, vs a monolithic HMXB en- 
semble. 

(3) Search for coronal X-ray emission from the most 
active stars in the SMC. 

(4) Constrain the incidence of LMXBs and CVs in the 
SMC. CVs should be present throughout the SMC which 
is dominated by an old population of red giants. Out- 
bursts (Novae and dwarf novae) are extremely infrequent 
and fainter in X-rays than HMXB. Thus they will be 
harder to find. No LMXBs or CVs are known in the 
SMC, and extrapolation from the total mass in stars pre- 
dicts only a handful of LMXBs , thus an observational 
constraint (or positive discoveries of LMXB) will provide 
important input to population synthesis models. 
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TABLE 1 
Chandra Observations 



Field 


ObsID 


Aimpoint 


Roll 


Start 


Exposure 






RA, Dec 


degrees 


UT 


xlO^ sec 


Deep Field 1 


7155 


00:53:34.50 -72:26:43.20 


161 


Apr 25 2006 05:32 


50 




7327 


00:53:34.50, 72:26:43.20 


161 


Apr 26 2006 15:14 


50 


Deep Field 2 


8479 


00:50:41.40 -73:16:10.30 


317 


Nov 21 2006 12:21 


45 




7156 


00:50:41.40 -73:16:10.30 


317 


Nov 22 2006 19:06 


39 




8481 


00:50:41.40 -73:16:10.30 


317 


Nov 23 2006 16:06 


16 
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TABLE 2 
Pulsars 



Name Source (CXO) srcid Period{s) PulsePraction Signif(%) Counts RA° Dec° r95%(") 



New Pulsars in Outburst 



SXP326" 


J005252.2-721715 


DF1_025 


326.79±0.16 


0.58± 0.07 


100.0 


9292 


13.217712 


-72.287619 


0.484 


SXP892'' 


J004929.7-731058 


DF2_019 


894.36±0.98 


0.65±0.15 


100.0 


1433 


12.374103 


-73.182874 


0.485 


New Pulsar Candidates 




J005331. 7-722240 


DF1_018 


131.114±0.025 


0.73± 0.26 


91.5 


386 


13.382157 


-72.378046 


0.373 




J005437.1-722637 


DF1_003 


10.4258±0.00016 


0.64±0.26 


90.4 


433 


13.654878 


-72.443838 


0.415 




J005337.8- 722409 


DF1_015 


7.83356±0. 00009 


0.78±0.33 


99.2 


160 


13.407884 


-72.402501 


0.353 


h 


J005448.9- 722544 


DF1_034 


358.03±0.19 


0.87±0.34 


90.8 


134 


13.704161 


-72.429125 


0.581 


b 


J005411. 5- 723513 


DF1_121 


665.32±0.65 


0.84±0.32 


94.6 


166 


13.54830 


-72.586961 


0.998 


b 


J005331. 8- 722646 


DF1.017 


594.53±0.52 


1.02±0.42 


97.2 


58 


13.382729 


-72.446133 


0.339 


a 


J005057.1-731008 


DF2_003 


98.774±0.012 


0.54±0.26 


71.6 


700 


12.737953 


-73.168905 


0.463 


a 


J004905.1-731411 


DF2_111 


8.9413±0.0001 


1.15±0.33 


93.3 


147 


12.271583 


-73.236528 


0.744 


a 


J005215.4- 731915 


DF2_022 


7.09496±0.00006 


0.615± 


61.8 


286 


13.064582 


-73.32096 


0.610 




J005059.0- 732055 


DF2_044 


9.3912±0.0001 


0.95±0.37 


96.2 


109 


12.746216 


-73.348768 


0.494 




J005042.2- 732137 


DF2_055 


13.2490±0.0002 


0.85±0.34 


98.8 


100 


12.675849 


-73.3604919 


0.553 




J004921. 0-731740 


DF2_077 


8.88999±0.0001 


0.10±0.40 


92.8 


68 


12.337766 


-73.294586 


0.720 




J005248.8-731730 


DF2_116 


23.5405±0.0007 


1.05±0.39 


91.1 


58 


13.203712 


-73.291884 


1.529 


c 


J005143.5-730408 


DF2.S 


6.30716±0.00005 




91.1 


153 


12.931325 


-73.068994 


2.132 


Known Pulsars Detected With Pulsations 


SXP138 ° 


J005323.8- 722715 


DF1_020 


138.925±0.028 


0.61±0.16 


100.0 


1232 


13.34958 


-72.4543259 


0.296 


SXP59 " 


J005456.1-722648 


DF1_032 


58.8336±0.0051 


0.82±0.21 


100.0 


398 


13.734005 


-72.446703 


0.504 


SXP756 " 


J004942.0- 732314 


DF2_017 


746.24±0.68 


0.52±0.10 


100.0 


4116 


12.42503 


-73.387348 


0.451 


SXP323 " 


J005044.6- 731605 


DF2.010 


317.26±0.12 


0.54±0.14 


100.0 


2311 


12.686208 


-73.268118 


0.289 


SXP172 " 


J005151. 9-731033 


DF2_001 


171.851±0.036 


0.39±0.17 


100.0 


3322 


12.96638 


-73.176053 


0.443 


SXP8.89 " 


J005153.3-723148 


DF1.S 


8.89917±0.00008 




99.3 


187 


12.972413 


-72.530032 


1.185 


Known Pulsars Detected Without Pulsations: 


SXP15.3 


,1005214.0-731918 


DF2_023 








147 


13.058339 


-73.321724 


0.711 


SXP7.78 


J005205.6- 722604 


DF1_153 








76 


13.023633 


-72.434479 


0.769 


SXP9.13 °- 


J004913.5-731138 


DF2_079 








1178 


12.306384 


-73.193905 


0.527 



Known Pulsars in the Field Not Detected: 



SXP4.78 DF-1 
SXP46.6 DF-1 
SXP82.4 DF-l^ 



"Quantiles indicate powcrlaw r=l spectrum with S/N>10 
''detected in binned X 16 lightcurve 
"^falls on ACIS-S detector 



TABLE 3 
Spectral Fits 



Source 


nH 


PL 


Spectral bins 


^red 


FX2-10 




10^1 cm 






erg/cm?/s 


57155B3-003 


1.57 ± 1.25 


1.81 ± 0.21 


26 


1.4 


6.668C-14 


57155B3-018 


2.39 ± 1.52 


1.84 ± 0.22 


23 


1.34 


6.290C-14 


57155B3-020 


2.30 ± 0.77 


0.76 ± 0.09 


71 


0.77 


5.241C-13 


57155B3-025 


3.44 it 0.23 


1.14 ± 0.03 


322 


1.25 


2.934e-12 


57155B3-032 


1.69 ± 1.54 


0.90 ± 0.21 


24 


1.39 


1.369e-13 


57155B3-042'' 


1.65 ± 0.94 


2.91 ± 0.30 


17 


0.24 


1.493e-14 


58479B3-001 


1.58 ± 0.29 


1.14 ± 0.05 


159 


1.02 


1.024e-12 


58479B3-003 


5.20 ± 1.48 


1.07 ± 0.1 


42 


0.86 


3.012e-13 


58479B3-010 


3.49 ± 0.59 


0.89 ± 0.07 


123 


0.98 


1.033e-12 


58479B3-017 


5.95 ± 0.49 


1.04 ± 0.05 


197 


0.98 


1.797e-12 


58479B3-019 


7.69 ± 1.01 


0.96 ± 0.10 


81 


0.96 


6.958e-13 


58479B3-022 


0.65 ± 1.50 


0.93 ± 0.20 


20 


1.21 


1.133C-13 


58479B3-074 


4.95 ± 2.14 


1.73 ± 0.24 


27 


1.46 


l.llle-13 


58479B3-079 


10.13 ± 1.33 


1.34 ± 0.12 


67 


1.27 


4.635e-13 



"This source is not a pulsar 
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Fig. 1. — Image of Deep Field 1 Chandra data, showing locations of known (boxes) and detected (circles) pulsars. North 
is up, East to the left. The ACIS-I field measures 16'xl6' on each side. The brightest source in the image is the new 
Pulsar SXP326. The many unidentified point sources are a mixture of quiescent X-ray binaries and background AGN. 



SMC Deep Fields 
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Fig. 2. — Image of Deep Field 2 Chandra data, showing locations of known (boxes) and detected (circles) pulsars. 
North is up, East to the left. The ACIS-I field measures 16'xl6' on each side. Note the two supernova remnants, one is 
the brightest object in the image, the other is much fainter and located to the upper right. Both SNRs have a Pulsar 
candidate within their boundaries. Pulsar SXP15.3 has a close (but resolved) neighbor which is itself a candidate pulsar. 
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Fig. 3.— SXP894. New Pulsar. P=894.36s 




Fig. 4.— SXP756. XS07156B0_014, p=746.25s 




Fig. 5.— SXP326. New Pulsar. XS07155B0_004, p=326.8s 



SMC Deep Fields 
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Fig. 7. — SXP172. Detected pulse period=171.8s The additional peak in the PDS at 235.8s is exactly half the Chandra 
ACIS Z-axis dither period of 707s. 




Fig. 8. — SXP138. Detected pulse period = 138.9s. The source was close to the edge of the ACIS detector, causing the 
additional peaks in the PDS which are artifacts of Chandra's dither pattern. 
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Fig. 9. — SXP98. New pulsar. The source J005057. 1-731008 had 700 net counts, and was examined for periodicities 
after it was found to exhibit an X-ray spectrum typical of SMC pulsars. 




Fig. 10.— SXP59 , p=58.8s 



31- 





Fig. 11.— SXP8.89. p=8.89s 



SMC Deep Fields 
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Fig. 13. — Candidate pulsar, with period 7.09s. The source CXO J005215.4-731915 had 286 net counts, and was 
examined for periodicities after it was found to exhibit an X-ray spectrum typical of SMC pulsars. 
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Fig. 14. — SMC Deep Field 1. Quantile Diagrams for cuts in Signal-to-Noise Ratio of >3, >5 & >10 (from top to bottom). 
Pulsars with periods were detected above 90% significance are indicated by black diamonds. Spectral model grids are 
shown for Power- law plus absorption (Blue) and thermal bremsstrahlung (Red). The confirmed pulsars cluster along 
the PL=1.0 contour and all are consistent with nH<10^^. In the upper left are several sources with small errors, these 
should be stars with coronal activity. In the bottom panel two grey triangle-points show the flaring source described 
in section [6l The upward triangle is the value during the flare, and the downward triangle is after the flare. The black 
point close to the flaring value is the result for the stacked data. 
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Fig. 15. — SMC Deep Field 2. Quantile Diagram. Pulsars whos' periods were detected above 90% significance are 
indicated by black diamonds. The black dots with error bars are plotted for every "good" point-source in the stacked 
100 ksec observation of Deep Field 2 . Spectral model grids are shown for Power-law plus absorption (Blue) and thermal 
bremsstrahlung (Red). The confirmed pulsars cluster along the PL=1.0 contour and all are consistent with nH<10^^. In 
the upper left are several sources with small errors, these should be stars with coronal activity, it will be very interesting 
to discover if any are in the SMC -as opposed to boring foreground galactic stars. 
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Fig. 17 — Power Spectrum of the flare in source CXO J005428.9- 723107. The input lightcurve is 7000sec of unbinned 
data containing the event. The PDS is plotted in log-log scaling to look for red-noise and QPOs. No identiflable features 
can be seen, indicating that the source is not a pulsar or magnetar. 




Fig. 18. — Spectrum of the flaring source CXO J005428. 9- 723107. The fltted model is a Raymond-Smith model (optically 
thin thermal plasma) with plasma temperature kT= 2.5 +/- 0.4 keV 



